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We consider a driven quantum particle in the strong friction regime described by tlie quantum 
Smoluchowski equation. We show that Crooks and Jarzynski like relations hold when the entropy 
production is properly generalized to take into account the non-Gibbsian character of the equilibrium 
distribution. We further propose a Josephson junction experiment that would allow to investigate 
nonequilibrium entropy fluctuations in overdamped quantum systems. 
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Thermodynamic processes at the nanoscale are gov- 
erned by both thermal and quantum fluctuations. It 
has lately been recognized that for classical nanosystems 
the second law of thermodynamics has to be generalized 
in order to include effects induced by thermal fluctua- 
tions. Thermal fluctuations are usually vanishingly small 
in macroscopic systems and are therefore neglected in 
the traditional formulation of thermodynamics pj . These 
generalizations of the second law take the form of fluc- 
tuation theorems that quantify the occurrence of neg- 
ative fluctuations of quantities like work, heat and en- 
tropy 0, 0] ■ A remarkable property of these new thermo- 
dynamic identities is their general validity arbitrary far 
from equilibrium. An important example of a fluctuation 
theorem is the one derived by Crooks, relating the prob- 
ability distributions of work, {W) and p^{W), along 
forward and reversed transformations of a system [3| , 



p'^{-W) = p'^iW) exp {-(3(W 



AF)) 
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Here AF is the free energy difference between final and 
initial states of the system. Equation ([T]) indicates that 
large negative work fluctuations are exponentially sup- 
pressed and hence not observable in macroscopic systems. 
In its integrated form, the Crooks relation reduces to an 
equality previously obtained by Jarzynski 



(exp i-P{W - AF))) = 1 



(2) 



connecting the equilibrium free energy difference AF to 
the nonequilibrium work W. In the above equation the 
average (...) is taken over the forward work distribution. 
It is important to realize that Eqs. ^ and ^ only 
apply for systems that are initially in an equilibrium 
Gibbs state. Extensions of these expressions for differ- 
ent nonequilibrium initial distributions have been intro- 
duced by Hatano and Sasa and by Seifert 0] . Fluctu- 
ation theorems have been investigated experimentally in 
various nonequilibrium situations [1, 0, [13, [HI, El, [3], 
the canonical example consisting of a highly damped 
Brownian particle in a driven potential. Due to the ex- 
perimental and theoretical importance of the strongly 
damped regime, the overdamped Langevin equation and 
the equivalent Smoluchowski equation have become the 



tool of choice for the analysis of classical fluctuation the- 
orems. 

In this paper, we derive quantum extensions of the clas- 
sical Crooks and Jarzynski relations, Eqs. ([T]) and ([2]), 
in the strong friction regime. Previous studies on iso- 
lated or weakly damped quantum systems can be found 
in Refs. |lU]. In the following, we use the quantum gen- 
eralization of the Smoluchowski equation to treat both 
thermal and quantum fluctuations 3, 121- Besides 
the known classical expression for the entropy produc- 
tion, we identify an additional contribution induced by 
quantum fluctuations. We further propose a Josephson 
junction experiment that would allow to test our predic- 
tions. 

Quantum Smoluchowski equation. In the limit of high 
friction, the off-diagonal matrix elements {x\p{t)\x') of 
the system density operator in the position representa- 
tion are strongly suppressed over a time scale of the order 
of 1/7, where 7 is the friction coefficient. As a result, a 
coarse-grained description of the dynamics of the system 
in terms of the diagonal part of the position distribu- 
tion alone, p{x,t) = {x\p{t)\x), becomes possible. In this 
semiclassical picture, the quantum system follows a clas- 
sical trajectory and quantum effects manifest themselves 
through quantum fiuctuations that act in addition to the 
thermal fluctuations induced by the heat bath. A notable 
advantage of this description is that the usual classical 
definitions of work and heat are valid, in contrast to the 
full quantum regime The quantum Smoluchowski 

equation can be written as [H, 3, 13] i 
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where V'{x) is the derivative of the external potential 
with respect to position and /3 the inverse temperature. 
The effective diffusion coefficient De{x) is given by, 

D,{x) = l/[l~X(3V"{x)] , (4) 

with the parameter 

A = (h/n-fm) [c + 277 /-f hp + * [27: /jh(3)] (5) 

measuring the magnitude of the quantum fluctuations. 
Here m denotes the mass of the system, c = —0.577... is 
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the negative Euler constant and ^' the digamma function. 
It should be noted that quantum corrections depend ex- 
phcitly on the position of the system through the curva- 
ture of the potential V"{x). When A = 0, Eq. ^ reduces 
to the classical Smoluchowski equation with constant dif- 
fusion coefhcient. The stationary solution of Eq. ([3]) is, 

Ps{x) = A^exp {-f3V{x) + Xl3^V'{xf/2) [1 - \pV"{x)] , 

(6) 

where N is the normalization constant. The above 
equilibrium expression is in general non-Gibbsian when 
A ^ 0. This observation indicates that the Crooks and 
Jarzynski relations cannot hold unmodified in the quan- 
tum regime and that proper generalizations arc needed. 

The quantum Smoluchowski equation ([3]) with the 
effective diffusion coefficient ^ is valid in the 
range of parameters, 7/^0 » ?i/3, 1/7, % 1//3 and 
\Xf]V"{x)\ < 1, where cuq is a characteristic frequency of 
the system [15|, |l6|, |l7|. In the present analysis, we as- 



sume that the quantum Smoluchowski equation is valid 
for time-dependent problems when the potential V{x, at) 
is driven by some external parameter at = a{t), that is, 
the non-diagonal elements of the density operator of the 
driven system remain negligibly small at all times. 

Quantum fluctuation theorems. We now derive exten- 
sions of the Crooks and Jarzynski relations, Eqs. ^ and 
(l2|), by using a path integral representation of the so- 
lution of the quantum Smoluchowski equation following 
Ref. 19]. For the sake of generality, we consider a generic 
driven Fokker-Planck equation, with position-dependent 
drift and diffusion coefficients Di and D2, of the form. 



di 



p{x,a,t) = Lap{x,a,t) , 



where the linear operator is given by, 
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The quantum Smoluchowski equation ^ corresponds 
to the particular choice Di{x^a) — --V'{x,a)/^m and 
D2{x,a) = 1/[1 — \(}V" {x,a)\"fm,f3. For any fixed value 
of the driving parameter a, we write the stationary solu- 
tion of Eq. ([7]) as, 



Ps [x, a) = cxp {-LP (x, a)) 



(9) 



where the function (p{x^ a) is explicitly given by 

ip{x,a) -- 



dy D2iy^^)-Di^iy^^) (10) 



D2 {y,a) 



We denote hy X = {x}^^^ a trajectory of the system 
that starts dX t = — r and ends at t = +t. We fur- 
ther define a forward process af , in which the driving 
parameter is varied from an initial value aZ-r = ao to 
a final value a^^ = c^i, as well as its time reversed pro- 
cess, a^ = a^^ . The conditional probability of observing 



a trajectory starting at x-t for the forward process can 
then be written with ut = {xt, af) as [20| . 



[X\x.r] = 



1 



D2 (Ur) 



exp 



- dtS [xt,xt,af) 



(11) 

with a similar expression for the reversed process. The 
prefactor in Eq. (jlip . which depends on the final value u-r, 
is essential in nonlinear systems to guaranty that thermo- 
dynamic potentials are independent of the state represen- 
tation [20| . In Eq. (fTTIl the generalized Onsager-Machlup 
function S {xt, if, a) is taken to be of the form 21[, 

o, ■ N [At ~ Di{xt,a)f 

S[xt,xt,a) = — — r . (12) 

AD2[xt,a) 

By assuming that the system is initially in an equilib- 
rium state given by the solution © of the Fokker-Planck 
equation, we obtain that the net probability of observing 
the trajectory X for the forward process is 



P^[X]^p, {u^r)P^ [X\X-r] ■ 



(13) 



In complete analogy, we find that the corresponding un- 
conditional probability for the reversed process reads 

(14) 

where we have introduced the time-reversed trajectory, 
X^ — {a;|}l^ with x\ = X-t- We next compare the prob- 
ability of having the trajectory X during the forward 
process with that of having the trajectory X^ during the 
reversed process. We have 



X^xl. 



1 



L>2 [U^r] 
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exp 
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dtS'' {xt,xt,af) , (15) 



where we have defined the conjugate Onsager-Machlup 
function, S't (^xt,Xt,Oif) — S (^Xt,—xt,af-). The ratio of 
the conditional probabilities pT|) and ([TS]) is simply de- 
termined by the difference of S and 5^. Using the defi- 
nition we thus obtain, 

[X\X-r] 



pR 



xn 



D2 (u-r) 
D2 (Mr) 



exp 



dt — ; — rXt 



D2 (ut) 



The ratio of the forward and reversed probabilities, 
Eqs. ((131) and (HH), follows directly as. 



P^[X] p,{u^r)P^[X\x^r] 
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where 



(18) 



By using the exphcit expression pO|) of the stationary 
solution ip{x, a), we finally arrive at 



= cxp 



[X] 

X exp / dt ^ ) [ it 



dt — at 

or 
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We are now in the position to derive generahzed fluctu- 
ation theorems for stochastic processes described by the 
generic Fokker-Planck equation ([7]). We begin by defin- 
ing the generalized entropy production E as, 



2^ — I at — — af 

oa 



at — — ; — r Xt -I- In ' 



D2 {ut) 



D2 {Ur) 



(20) 

The first term in Eq. (|20p is similar to the entropy pro- 
duction introduced by Hatano and Sasa for systems ini- 
tially in a nonequilibrium steady state [6|] . In the present 
situation, however, it corresponds to a non-Gibbsian 
equilibrium state. The second and third term, on the 
other hand, find their origin in the position dependence of 



the diffusion coefficient and vanish for constant D2. We 
note in addition that the entropy production as defined in 
Eq. ^ is odd under time-reversal, [Xt] === -E^ [X]. 
The distribution of the entropy production, p^(E), for an 
ensemble of realizations of forward processes can then be 
defined as, 

(E) = JvXP" [X] ^ (E - E^ [X]) ^ 
= exp (S) J VX^ [X^] (5 (E -t- E^ [X^] ) (21) 

where we have used Eq. in the last line. Here / VX 
denotes the product of integrals over all possible paths 
X. The continuous integral in Eq. pip is interpreted 
as the limit of a discrete sum in the Stratonovich sense. 
Equation ((2T|) can be recast in the form of a generalized 
Crooks relation for the entropy production, 



p^(-E)=p^ (E)exp(-E) 



(22) 



By moreover integrating Eq. (|22|) over E, we obtain an 
extended version of the Jarzynski equality. 



(exp(-E)) = l 



(23) 



For the case of the quantum Smoluchowski equation ([3]) , 
the generalized entropy production (|20p reads. 



dt' 



1- pXV"{xuaP) 



it -f In 



• I- P\V" {xr,ai) ' 
1~PXV" [x-r^ao], 



(24) 



Expression (|24p for the entropy production, together with 
the fiuctuation theorems ((22)) and (p3|) . constitutes our 
main result. Combined, they represent the quantum 
generalizations of the Crooks and Jarzynski equalities, 
Eqs. dl]) and ([2|), in the limit of strong damping. The 
first term in Eq. ([^H) is analogous to the classical en- 
tropy produced when the system is driven out of equilib- 
rium when at is varied in time. By contrast, the second 
and third terms are induced by quantum fluctuations and 
vanish in the classical limit A = 0. In the latter case, 
<p{u) = P{V{u) — F{a)) and the entropy takes the fa- 
miliar form, Yj = P j dtdV{ut)/daaf — (3AF. We note 
that the total entropy production is zero for undriven 
systems, in agreement with detailed balance (the non- 
trivial prefactor in Eq. (jlip is paramount to obtain this 
result). The inequality (E) > implied by Eq. ^ is 
often interpreted as an expression of the second law. 

Experimental verification in Josephson junctions. No 
experimental investigation of quantum fiuctuation theo- 



rems has been performed so far. A scheme to study the 
Crooks and Jarzynski relations in isolated and weakly 
damped quantum systems using modulated ion traps has 
recently been put forward in Ref. [23 |. Here we propose 
to test the predictions for the quantum fluctuation the- 
orems, (1^^ and in the strong damping limit using 
Josephson junctions [23|. The Josephson relations for the 
current Is {t) across the junction and the phase difference 
(j)(t) between left and right superconductors are 

Is = Ic sin ((/>) and <j)^2e/h U{t) , (25) 

where U{t) is the voltage drop across the junction. The 
maximal current Ic is given by Ic = {2e/h)Ej^ where 
Ej is the coupling energy (Josephson energy). An ex- 
ternally shunted Josephson junction can be described 
via an equivalent circuit consisting of an ideal junc- 
tion, a capacitance C and a resistance R (Resistively 
Shunted Junction (RSJ) model) [2^. In this model, the 
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TABLE I: Typical parameter values for Josephson junctions. 





T[K] 


C[pF] 


R[n] 


Ic[mA] 




A[10-^] 


e 


Classical 


4.2 


0.1 


0.37 


1.0 


0.041 


-0.0067 


-0.038 


Quantum 


4.2 


0.1 


1.0 


1.0 


0.30 


-1.1 


-6.3 



Josephson junction is interpreted as describing the dif- 
fusive motion of a particle with position (j){t) and mass 
m — {h/2e)'^C, the friction coefficient being given by 
7 = 1/RC. An important quantity in the RSJ model is 
the dimensionless capacitance (Stewart-McCumber pa- 
rameter), /3c = {2tt /^q)IcR^C, where $o = h/2e is the 
magnetic flux quantum. In the overdamped regime, 
/3c < 1, the dynamics of the Josephson phase can be de- 
scribed by the quantum Smoluchowski equation ([3]) with 
the potential V{(j)) = -Ejcos{4>) - Ejcj) '^]. The en- 
ergy Ej = [h/2e)I is determined by the bias current / 
and the effective diffusion coefficient Q reads, 



De{4,) = l/[l-e cos (,/))] 



(26) 



The constant Q = X(3Ej is the crucial parameter which 
governs the magnitude of quantum effects in a Joseph- 
son junction. It is directly proportional to the quantum 
parameter A, Eq. ([5]), which in the context of the RSJ 
model can be reexpressed as. 



X = 2p[c + 2Tr^ p/(3Ec + * {/3Ec/2tt^p)] 



(27) 



where Ec — 2e^/C is the charging energy, p ~ R/Rq 
the dimensionless resistance and Rq = h/Ae^ the resis- 
tance quantum. In Tab. U we list typical parameter val- 
ues for niobium-based Josephson junctions. By varying 
the value of the resistance, both the classical, 6^1, 
and the quantum regime, 6^1, can be realized. 

The nonequilibrium entropy production (j24p can be ex- 
perimentally determined in a Josephson junction by ap- 
plying the following measurement procedure. The phase 
4> can be directly deduced from a measurement of the 
Josephson current once the current -phase relation of the 
junction has been determined [25|. The system is then 
first prepared in a given initial state and let to relax to its 
stationary state (O. After the latter has been attained, 
the Josephson potential V{(j)) is modified according to a 
specific driving protocol at with the help of an external 
magnetic field. The entropy production E during such a 
protocol (corresponding to either a forward or reversed 
transformation) can be evaluated via Eq. ([24|) from the 
recorded values of the current. The distribution function 
of the entropy can eventually be reconstructed by repeat- 
ing the above measurement sequence and the validity of 
the quantum fiuctuation theorems ((22)) and ([23]) in the 
strong friction regime can be tested. 

Conclusion. We have analyzed quantum fluctuation 
theorems in the strong coupling limit with the help of 
the quantum Smoluchowski equation. We have shown 



that the Crooks and Jarzynski relations remain valid in 
this regime when the entropy production is properly gen- 
eralized to take into account the non-Gibbsian property 
of the equilibrium state. We have moreover proposed 
a Josephson junction experiment that would enable to 
study quantum nonequilibrium entropy fluctuations in 
overdamped systems. 
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